Quasi-bound states of Schrodinger and Dirac electrons in magnetic quantum dot. 



o 
o 

(N 



M. Ramezani Masir/'Q A. Matulis, 1, 2 'Q and F. M. Peeters 1 ' 3 '^ 

1 Departement Fysica, Universiteit Antwerpen 
Groenenborgerlaan 171, B-2020 Antwerpen, Belgium 
2 Semiconductor Physics Institute, Gostauto 11, LT-01108 Vilnius, Lithuania 
5 Departamento de Fisica, Universidade Federal do 
Ceard, Caixa Postal 6030, Campus do Pici, 60455-760 Fortaleza, Ceard, Brazil. 

(Dated: July 21, 2009) 

The properties of a two-dimensional electron are investigated in the presence of a circular step 
magnetic field profile. Both electrons with parabolic dispersion as well as Dirac electrons with linear 
dispersion are studied. We found that in such a magnetic quantum dot no electrons can be confined. 
Nevertheless close to the Landau levels quasi-bound states can exist with a rather long life time. 
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I. INTRODUCTION 

During the last five years, graphene (a single layer of 
carbon atoms) has become a very active field of research 
in nanophysicsii 2 -. It is expected that this material will 
serve as a base for new electronic and opto-electric de- 
vices. One of the most challenging task is to learn how 
to control the electron behavior using electric fields in 
this two-dimensional (2D) layer. This task is made com- 
plicated by the so called Klein effect according to which 
Dirac electrons in graphene can tunnel through any elec- 
tric barrier—. As a consequence in electrically created 
quantum dots there are no bound states but only quasi- 
bound states, or so called resonances^^ which, however, 
under certain conditions can have a long life time. 

An alternative approach to control the motion of elec- 
trons is to use non homogeneous magnetic fields which 
can be created e.g. through the deposit of nonostruc- 
tured ferromagnets. 7,8,9 . Recently it was shown that 
nonhomogeneous magnetic structures are able to con- 
fine Dirac electrons in graphen o 10 ' 11 ! 12 : 13 ! 14 : 15 ! 16 ! 17 . How- 
ever, up to now semi-infinite (homogeneous in one direc- 
tion) structures were considered, what makes the anal- 
ysis more simple because the problem is reduced into a 
one-dimensional one. 

In this paper we consider a finite size magnetic struc- 
ture where the magnetic field is nonzero only in a finite 
region of space. Namely, we consider a model homoge- 
neous magnetic field that is non-zero in a circle that we 
call the magnetic dot. This situation is the inverse of 
the one considered in Ref. 11 where a magnetic anti-dot 
was considered, as in Ref. 9 for the case of normal elec- 
trons, where the magnetic field is zero in a circular region 
and non-zero outside this region. Such a model system 
can be realized by having a magnetic vortex piercing the 
graphene layer or by overlaying graphene with type I su- 
perconductor with a circular hole placed in the perpen- 
dicular magnetic field. In order to reveal the peculiarities 
of the behavior of Dirac electrons in such magnetic dot 
we compare the result with those for standard electrons 
with parabolic dispersion law. 



We show that it is impossible to confine 2D electrons 
in a magnetic dot in contrast to semi-infinite magnetic 
structures neither in the case of graphene nor in the case 
of the standard electron, and consequently, all Landau 
levels convert themselves into unbound states. Never- 
theless, long living quasi-bound states can be present. 
We studied them using the local density of states tech- 
nique applied previously for the investigation of electri- 
cally confined electrons^. 

The paper is organized as follows. In section II the 
model of a magnetic dot for a standard electron is con- 
sidered. The problem is formulated in subsection A, the 
local density of states technique is presented in subsection 
B, the results are discussed in subsection C, while in sub- 
section D the complex energy eigenvalues of the problem 
are described. In the corresponding subsections of sec- 
tion III the problem of the Dirac electron in a magnetic 
dot is presented. Our conclusions are given in section IV. 



II. ELECTRON WITH PARABOLIC ENERGY 
DISPERSION 

We assume that a homogeneous magnetic field Bo is 
present in a circular area of radius tq , while there is no 
magnetic field outside it, namely, Brj(r) = e z Bo®(ro — r). 
The behavior of the electron is described by the station- 
ary Schrodinger equation 



{H- £}*(r) = 0, 
with the Hamiltonian 



(1) 



(2) 



Because of the cylindric symmetry of the problem we 
choose the symmetric gauge for the vector potential 
defining its single azimuthal component as 



1 



r < r ; 



Aip{r) — 2 \ rg/r, r < r. 



(3) 



This azimuthal component is shown in Fig. Q] together 
with the magnetic field profile. In order to simplify the 
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FIG. 1: (Color online) Azimuthal vector potential component 
A v (blue solid curve), and perpendicular magnetic field B z 
(red dashed curve) as functions of the radial coordinate. 



notations we use dimensionless variables, based on the 
magnetic field strength value B$. Thus, the magnetic 
field B(r) is measured in Bo units, all distances are mea- 
sured in the unit of magnetic length Ib — y^ch/eBo, 
energy and potential in fko c {oj c = eBo/mc), and vector 
potential in BqIb units. In the case of electron moving 
at a GaAs/AlGaAs interface (m* = 0.067) and magnetic 
field of 1 T the unit of length is Ib = 250 nm, and the 
energy unit is 20 meV. 



A. Solution of eigenvalue problem 

The Schrodinger equation |T]) in cylindric coordinates 
reads 



ld_ d 



1 d 2 



iA v d 



r dr dr r 2 d(p 2 r dip 
Substituting the wave function 

* = *(r,p) = e im 
we arrive at the radial equations 

2 



A, 



2E \^ = 0. 



(4) 
(5) 



1 d d 

r dr dr 



(™ + £) +2E^ I {r)=Q, (6a) 
Id d (m + rl/2) 2 1 

~^ r ^~- § l - l + 2E\^ ii {t) = q, (6b) 

r dr dr r z 



which have to be solved inside the dot (region I) and 
outside it (region //). The boundary conditions (the 
continuity of the wave function and its radial derivative) 
have to be satisfied at the dot border (r = ro). 

The regular solution inside the dot can be expressed via 
the confluent hypergeometric function (Kummer function 
M(a\c\z)): 



^j(r) =Af(r) =Ar^e~ r2 ^ 

x M ((|m| +m)/2 + 1/2 - B||m| + l|r 2 /2) , 



(7) 



while the solution outside it is composed of two Bessel 
functions 



where k = y/2E is the momentum of the free electron 
(measured in Ig 1 units), and v — m + Tq/2. Note both 
functions ( J„ and Y v ) suit us, as they vanish in the limit 
r — > oo. 

Thus, we have three constants A, B, and C . They can 
not be defined from the above mentioned two boundary 
conditions. That is why we have to conclude that there 
are no bound states, and consequently, a magnetic field 
in a finite region of the 2D plane can not confine the 
electron. However, quasi-bound states can be expected 
when the electron energy in the dot is close to the Landau 
levels with energy 



E,, 



1 



(9) 



(here n = 0,1,--- and m = 0,±1,---) defined in the 
case of homogeneous magnetic field. Confirmation of this 
statement follows from Fig. [2 where the electron wave 
functions for two different energies are shown. 




^ u {r) = BJ v {kr)+CY v {kr), 



(8) 



FIG. 2: (Color online) The wave functions for m = 0, ro = 5: 
E — 2.5 - red solid curve, and E = 2.8 - dashed blue curve. 
The magnetic dot region is indicated by shadowed rectangle. 



We see that in the case of E = 2.5 (red solid curve) 
which corresponds to the Landau level with n = 2 and 
to = the wave function is large in the dot region (shown 
in Fig.[2]by shadowed yellow rectangle), while for the case 
of energy E = 2.8 , which does not coincide with any 
Landau level energy, it does not have any appreciable 
large value inside the dot, and actually does not differ 
much from the wave function for a free electron calculated 
in cylindric coordinates. 



B. Local density of states 

Next we will look for possible long living quasi- 
stationary states in the magnetic dot. In principle such 
quasi-bound (or quasi-stationary) states have to be de- 
scribed by the solution of the time dependent Schrodinger 
equation which is much more complicated as compared 
with the standard eigenvalue problem. There are, how- 
ever, several alternative approaches which enable us to 



3 



investigate properties of quasi-bound states by station- 
ary means. We follow the method presented in detail 
in Ref. and calculate the local density of states. The 
basic idea is to confine the electron in a large region of 
finite radius R, where its wave function obeys the zero 
boundary condition at the border (r = R) and treat the 
problem as a stationary one. A measurement that probes 
quantum dot properties, say, measuring of the tunnelling 
current directed perpendicular to the dot with STM, or 
power absorption in near field infrared spectroscopy, has 
to depend on the averaged value of the electron wave 
function in the dot. Therefore we introduce the integral 

/>oo 

1(E) =2n rdrF(r)\^(r)\ 2 , (10) 
Jo 

which depends on the electron wave function, and actu- 
ally is proportional to the so called local density of states. 
The aperture function F(r) characterizes the interaction 
of the electron with the measuring probe. 

This integral is sensitive to the probability to find the 
electron in the dot, and in the case of a quasi-bound state 
it will exhibit a peak corresponding to the energy of this 
state. The width of the peak is related to the inverse of 
the life time of this quasi-stationary state. 

For the sake of determinacy we use the aperture func- 
tion of a gaussian: 

F(r) = br 2 e~ br \ 6 = rjjT 2 In 10. (11) 

which corresponds to the probability to find the electron 
in the dot area Trr 2 ,. In the case of larger b value instead 
of the local density of states we obtain the squared wave 
function value in the center of the dot, while in the case 
of smaller b value the peculiarities of the dot are washed 
out. 

The solution of the Schrodinger equation ^ given by 
Eqs. (J7J) and © has to satisfy the following boundary 
conditions: 

Mro) = Mr ), (12a) 
4>i,r(r ) = ipn, r {ro), (12b) 
Tpn{R) = 0, (12c) 

which converts our problem into an eigenvalue problem. 
Here and further the subscript ,. means the derivative 
over r. 

At the end, we are interested in the limiting case R — ► 
co. Therefore, in the last of Eqs. (fT2"]) we replace the 
Bessel functions by their asymptotic, namely, we have 

B cos(kR — ip m ) + C sin(fci? — <p m ) — 0, 

ip m = 7T {to + (rg + l)/2} /2. ""' 

instead of Eq. (|12cp . 

Postponing till later the proper wave function normal- 
ization we assume that B = cos<l> and C = sin<i> and 
rewrite the above equation as 

cos{kR - tp m - $) = 0. (14) 



This equation shows that the eigenvalues of the consid- 
ered problem are approximately separated by Afc = tt/R, 
and reduce to a continuum spectrum in the limit R — > oo. 
Constructing some averaged description which is valid 
when calculating the local density of states, we replace 
Eq. (|12c[) by the following one: 

B 2 + C 2 = l. (15) 

Now solving it together with Eqs. (|12al b) we obtain three 
constants: 

A = —W, B = WQ, C = -WP, (16) 

TTTq 

with 

P = Jufr-Ju,rS, (17a) 

Q = Y u f r -Y Vtr f, (17b) 

W = (P 2 + Q 2 y 1/2 . (17c) 
The obtained constants enable us to calculate the integral 

CJJ. 

In order to convert the above integral into the local 
density of states we have to multiply it by two additional 
constants. One of them is the wave function normaliza- 
tion factor N, which can be estimated calculating the 
integral of the squared wave function in the limit of large 
radius R. The replacement of the Bessel functions by 
their asymptotic immediately leads to N = k/2R, The 
second one is a consequence of the replacement of the 
summation over the discrete eigenvalues by the integra- 
tion over energy, which is given by the factor R/irk. To- 
gether they give 1 /2ir, which results into definition of the 
local density of states 

p(E) = ±1(E). (18) 

Z7T 

C. Numerical results 

We solved numerically Eqs. (jTTJ) . Inserting the ob- 
tained results into (|16p. and later in Eqs. {7J and ([8]) we 
obtained the wave function what enabled us to calculate 
the integral (fTQ| . and finally the local density of states 

A typical result for the local density of states as a func- 
tion of electron energy is shown in Fig. [3J We clearly see 
peaks close to the energies of the Landau levels © cal- 
culated for the case of a homogeneous magnetic field. 
These peaks are broadened indicating that they are not 
really bound states in the magnetic dot. The broadening 
is larger for higher energy peaks. 

The next thing which also is seen in Fig. [3] is a decreas- 
ing background with energy. This background is due to 
the states of the free electron in the absence of the mag- 
netic dot. To justify this statement we made the same 
averaging (over circle of radius ro) with the same gaus- 
sian aperture function pip of the radial component of the 
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FIG. 3: (Color online) The local density of states for m — 
and ro = 3 shown by the solid red curve. The same density 
calculated for a free electron according to Eq. (|20|l is shown 
by the blue dashed curve. 



free electron wave function (when there is no magnetic 
dot). This function reads 

4>fr cc (r) = J m (kr) (19) 

and is valid in the whole 2D plane. Inserting this function 
into integral (fTTJj) and later in Eq. (|18p and using Tables 
of integrals^ we obtain the local density of states for a 
free electron 

Pf[ee (E) = r l e - E / b I m {E/b), (20) 

where I m {x) stands for the modified Bessel function of 
the first kind. This local density of free electron in the 
case of m — is shown in the same Fig. [3] by the blue 
dashed curve. Comparing these two curves we clearly see 
how increasing the electron energy we reduce the influ- 
ence of the magnetic dot on the electron behavior, and 
the local density of states converts itself gradually into 
the free electron one. 

We fitted the peaks in the density of states by 
Lorentzian functions a n ~f n /{(E — E n ) 2 + 7^} defining 
three parameters for any of them: the position E n , its 
broadening 7„, and the amplitude a n . Two of them (the 
position and broadening) are shown in Figs. 0] and [5] for 
different orbital momenta as functions of the radius of 
the dot r . The position of the quasi-bound states E n 
are shown by the red solid curves while the broadening 
of the peaks is indicated by the shadowed areas limited 
by the E n ± 7„ curves. 

Notice that the levels to the right of the green dotted 
curve are extremely narrow and their position coincides 
with the Landau levels © shown by the blue dashed hor- 
izontal lines. In fact this means that almost all electron 
wave function is located in the magnetic dot (using the 
classical description language we may say that the elec- 
tron rotates along the Larmor circle inside the dot) and it 
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FIG. 4: (Color online) Quasi-bound states with orbital mo- 
mentum m = 0. The energies of these states are given by 
red solid curves and the widths (i. e. the inverse of the life 
time) by shadowed regions. The Landau levels are indicated 
by blue dashed lines. 
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FIG. 5: (Color online) The same as Fig.f4]but now for m = 1. 



does not touch the border of the magnetic dot. When the 
dot radius r becomes smaller the Larmor circle touches 
the dot border and tunnelling of the electron outside the 
dot starts which broadens the level. The partial pene- 
tration of the wave function outside the dot leads to a 
lowering of the quasi-bound state energy as well. The 
raising of this energy for small ro values is caused by 
the large asymmetry of the peak where actually the ap- 
proximate replacement of the peak by a Lorentzian type 
function is no longer valid. This picture is more or less 
the same for all positive m values (compare Figs. 0] and 
[5J. The difference is that for larger m values the levels 
start at higher energies what is in agreement with the 
expression for Landau levels ([9|). 
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The picture for negative m values is different as shown 
in Fig. [51 All of them belong to the same Landau level 
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m = -1 




FIG. 6: (Color online) The lowest quasi-bound state with 
n = and different negative m values. The vertical dotted 
green lines are the analog of the dotted curve in Figs. [4] and [5] 
separating the weakly broadened states from those with small 
lifetime. 



energy which is an expression of the degeneracy of the 
Landau level. We see that with decreasing radius of the 
dot ro the levels with different m disappear step by step, 
the ones with smaller absolute m values disappear later. 
This is in agreement with the fact that the larger the |m| 
value the larger the radius of the electron trajectory, and 
the electron wave function is closer to the dot edge. 

The increase of the peak broadening at small ro values 
is so steep that it is worth to divide all the peaks into two 
classes as shown in the above figures by the green dotted 
curves. The levels on the left side of these curves belong 
to essentially broadened quasi-bound states, while those 
on the right side from the experimental point of view can 
hardly be distinguished from the real bound states. 

One can rudely estimate the position of this dividing 
curve comparing the approximate dimensions of the elec- 
tron wave function calculated in the case of a homoge- 
neous magnetic field (which actually coincides with func- 
tion ([7]0 with the magnetic dot radius ro. A more ac- 
curate estimation can be obtained solving the stationary 
Schrodinger equation for complex energy eigenvalues as 
was described in Ref. 0. We draw these dividing curves 
in Figs. HUH] using this technique which is sketched in the 
following section. 



D. Complex energy technique 

According to Ref. the lifetime of the quasi-bound 
state, or the trapping time of the electron in the quan- 
tum dot can be estimated solving the time-independent 
Schrodinger equation and applying boundary conditions 
of the outgoing waves at the sharp dot border. 

In our case this technique reduces to connecting the 
wave function ([7]) defined in the dot with the outgoing 



electron wave function outside it, which is given by the 
first kind Hankcl function 



^out(r) = H<p(kr) = J u (kr) + iY v {kr) (21) 



with 



v = m + Tq/2. 



(22) 



Applying the boundary conditions for these wave func- 
tions and their derivatives we obtain the following equa- 
tion: 



f r (r )H(kr )-f(r)H r (kr )=0, 



(23) 



where for the sake of simplicity we omitted the indexes 
of the Hankel functions. The indexes which are only left 
indicates the derivative over the coordinate r. This equa- 
tion has to be solved for a complex energy (or complex 
k), the imaginary part of the energy gives the inverse of 
the lifetime. 

For finding the dotted curves in Figs. HUH] separating 
the quasi-bound and nearly bound states it is enough 
to solve the above equation by means of a perturbation 
expansion in terms of the momentum difference Afc = 
k — ko where ko = y / 2E n , m with energy E n>m of the 
unperturbed Landau level. Limiting ourselves to first 
order in Afc we arrive at the following expression: 



Afc 



H r f — Hf r 



Hf r ,k + Hkfr ~ fH r , k - H r f k 



(24) 



All functions and their derivatives over r and fc have to 
be calculated at r = r$ and fc = ko. 

Now introducing the energy deviation from the Landau 
level energy 



AE = E- En 



fcnAfc, 



(25) 



taking its imaginary part and equating it to 10 -2 (it is 
expected that a smaller broadening can hardly be re- 
vealed experimentally) we obtained the points connected 
by the green dotted curve in Figs. [4][6] separating the 
quasi-bound states from nearly bound states. 



III. DIRAC ELECTRON IN GRAPHENE 

Now we repeat the above calculation for the mag- 
netic dot applying it to the case of a Dirac electron in 
graphene, where the low-energy quasi-particles (electrons 
and holes) are described by the following dimensionless 
Dirac-like Hamiltonian: 



H = a- (-iV + A) . 



(26) 



Here, <r — {a x ,a y } stands for the 2x2 Pauli matrices. 
The units are based on the magnetic field strength Bo and 
they are the same as in previous section, except the unit 
of energy which now is vf^/Ib with the Fermi velocity 
vp = 10 8 cm s _1 . In the case of a 1 T magnetic field this 
energy unit is 2.6 meV. The vector potential is given by 
Eq. ©. 
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A. Solution of eigenvalue problem 

The approach is based on the same stationary 
Schrodinger equation {1} but now with the matrix Hamil- 
tonian (j2l))) , which results into a set of two differential 
equations. Assuming the wave function of the following 
form: 



a[r) 



ie lv b(r) 



(27) 



we arrive at a set of two equations for the radial wave 
function components 



4- + A(r) + ^±1 \b = Ea, (28a) 
dr r 

fl IT) 

— - AM S> a = Eb, (28b) 

dr r 1 



which has to be solved in the two regions (I in the dot, 
and II — outside it). We require the continuity of the 
obtained components at the dot border ro 



ai(ro) = a H (r Q ), bi(r ) = b n (r ). 



(29) 



Instead of solving these first order differential equa- 
tions it is more convenient to convert them into second 
order differential equations for a single component, say 
for component b 



1 d d (m+ l) 2 



r- 

r dr dr 

1 d d 
r dr dr 



r 

T 



[E 2 



E 2 



(m+ 1 + r§/2) 2 



3/7 



0, (30a) 
: 0. (30b) 



In contrast to the electrical quantum dot case which 
was considered in Ref. [5| now the effective potential in 
Eqs. (f2"8"| is a continuous function at the dot border. For 
this reason the boundary conditions (f2"9"| are equivalent 
to 

bi(r ) = bu(r ), &i>(r ) = bu >r (ro)- (31) 

These boundary conditions are identical to those for the 
previous Schrodinger electron case (|T2|) . It enables us to 
use the full analogy with the previous case. Taking this 
analogy into account we have for the solution in the two 
regions 

bi(r) = Af(r) = ^rl m+1 l e - r2 / 4 M( ao |co|r 2 /2), (32a) 
b H (r) = BJ v (kr) + CY v {kr). (32b) 

where k = \E\, a = (|m + 1| + m + 1 - E 2 )/2, v = 
m + 1 + Tq/2, and cq — \m + 1| + 1. The expressions for 
the other wave function component a(r) follow directly 
from Eq. !|2"8"a|) : 

.4 



= _ r \m+l\ e -r 2 /i 

E 



ai{r) — —r ' '< 



|m + 1| + m + 1 



[dr r 
a H {r) = BJ„_i(fcr) + CY V - X {kr) 



M(a \c Q \r 2 /2), (33a) 



The wave function components obtained in the above 
way are illustrated in Fig. [7] for two different values of 
the energy. We see the same tendency. When the energy 



£ = 1 .932 




FIG. 7: (Color online) The wave function components: a - 
dashed blue curve, and b - red solid curve, for m = 0, dot 
radius ro = 5 and two energy values: E — 1.932 - upper plot, 
E = 2.470 - lower plot. 



is close to the Landau level energy of the Dirac electron 
in a homogeneous magnetic held, 



En 



±V2r 



1 



1. 



(34) 



(see the lower plot of Fig. [7] where the energy is close 
to the Landau level with m = 0, n = 2) we see a clear 
accumulation of the wave function components in the dot, 
what indicates a quasi-bound state. 



B. Local density of states 

Developing further the analogy of Eq. (|30]) with the 
considered previously case we calculated the local density 
of states using Eqs. (fT7) . (flU)) and (flT))) . Because we have 
now a wave function with two components, Eq. (|10[) is 
modified into 

poo 

1(E) = 2tt rdrF(r){\a(r)\ 2 + \b(r)\ 2 } . (35) 

Now the normalization factor is N = k/AR (due to the 
two wave function components), and the factor responsi- 
ble for the change of the summation over discrete eigen- 
values into an integral over the electron energy is R/ir. 
Thus, the local density of states in the case of a Dirac 
electron becomes 

P[E) = rdr f( r HWr)\ 2 + \b(r)\ 2 } ■ (36) 

In the case of free Dirac electrons (when there is no 
magnetic dot) the wave function components read 



Off 



J m (kr), 6frcc = J m +i(fcr), 



(37) 



(33b) 
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what leads to the following expression of the local density 
of states for a free electron: 

p bee (E) = \M e -EV2b { Im{E 2 /2b) + I m+1 {E 2 /2b)} . 

(38) 

C. Numerical results 

The typical local density of states calculated for m = 
and rp = 3 is shown in Fig. [H] for positive energies. Two 
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FIG. 8: (Color online) The local density of states for a Dirac 
electron in the magnetic dot for m = and ro = 3 shown by 
red solid curve. The dashed blue curve is the free electron 
density of states. 



differences with respect to standard electrons can clearly 
be noticed. First, in the case of the Dirac electron the 
spectrum is symmetric with respect to energy inversion 
(E — > — E) due to the equivalence of electrons and holes. 
Thus the plot in Fig. [8] has to be supplemented by the 
same curves for negative energies. Second, comparing the 
density of states for Dirac electron with the same curve 
for the Schrodinger one (see Fig. [3]) we see that there are 
more peaks. This can be explained by the more dense 
Landau level spectrum in the case of Dirac electrons (f34|) 
for the large quantum number values as compared with 
these for the previous case ©■ 

As before we fit the peaks by Lorentz type curves what 
leads to the broadened levels displayed in Figs. 151 and [ITU 

We see that the levels with m = — I start at lowest 
energies what is just the consequence of chosen definition 
of radial wave function components (|27|) . 

The green dotted curves divide the region of broadened 
quasi-bound states from the region where the states have 
a very small broadening. These curves were obtained in 
the same way as it was done in section III Dl for the case 
of Schrodinger electron, namely, applying the complex 
energy eigenvalue technique. In the case of the Dirac 
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FIG. 9: (Color online) Quasi-bound states with orbital mo- 
mentum m = for the Dirac electron in the magnetic dot. 
The energy of these states are given by red solid curves and 
its width (i. e. the inverse of the life-time) by the shadowed 
regions. The Landau levels are indicated by blue dashed lines. 
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FIG. 10: (Color online) The same as Fig. [4] but now f° r m = 
-1. 



electron it leads to the following imaginary energy part: 



H nu a r + H vr a — H v ^±b r — H v -i r b 

The above mentioned dotted green line corresponds to 
S = 3 • 10~ 3 . Note we chose it three times smaller than 
in Figs. H][6] what causes us to conclude that between 
the local density of states technique and complex energy 
eigenvalue method one can expect only qualitative agree- 
ment. 

There is one more interesting point — the zero energy 
state which is shown in Fig. [10] by thick red line along 
the x-axis. Its behavior differs essentially from all other 
states. That is why it needs some special attention which 
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is presented in the next subsection. 



D. Zero energy state 

Now we check whether the Dirac electron has a zero 
energy state in the magnetic dot. In this case instead of 
Eqs. (|28|) we have to solve the following two equations 
for the radial components of the electron wave function: 



— + AM + [b = 0, 

ar r 



Tr -A{r)--U = 0. 



(40a) 
(40b) 



These are uncoupled differential equations of the first or- 
der, and their solution can be found by a straightforward 
integration. The solution has the following asymptotic 
behavior: 



ln< 



(r) = [ dr{A{r) + ^ 



mlnr + r 2 /4, r — ► 0; 
(m + TQ/2)lnr, r — ► oo, 



(41) 



and 



lnfc(r) = - 



/*{ 



A(r) 



TO + 1 



r 

(to + 1) In r - r 2 /4, r — > 0; 
(to + 1 + Tq/2) lnr, r — » oo, 



. r m exp(r 2 /4), r — > 0; 



(42) 



and 



6(r) 



r- m - 1 exp(-r 2 /4), r -► 0; 

r -(m+l+r„/2) r ^ OO. 



(43) 



(44) 



In order to have the wave function with finite norm 
two boundary conditions have to be satisfied. First, the 
function should behave like r a (a 0) when r — > 0, and 
second, it should behave like r~ a (a ^ —1) when r — ► oo. 

For the a component the above conditions reduce to 
the requirements m and m + r 2 /2 ^ — 1, which can 
not be satisfied simultaneously. Consequently, we have 
to assume that a = 0. 

In the case of component b the conditions read 



r 2 /2 ^ m < -1, 



(45) 



from which it follows that if Tq/2 > 1 there are always 
some negative to values for which a zero energy state 
exists. When the radius of the dot decreases this interval 



becomes smaller, and the zero energy states vanish one 
by one. Finally, at r$/2 < 1 all of them disappear. 

Such essential difference between the bound zero en- 
ergy level and all other quasi-bound levels is caused by 
the fact, that the wave function of the state with zero 
energy is real. Consequently, the electron in this state 
has no velocity, and as a result there is no tunnelling of 
this electron outside the dot. Unfortunately, the absence 
of any non zero electron velocity makes it impossible to 
reveal this state in transport measurements, but maybe 
it can reveal itself through the statistic properties of the 
magnetic dot. 



IV. CONCLUSIONS 

We considered the eigenvalue problem of a model quan- 
tum magnetic dot where the homogeneous magnetic field 
perpendicular to the 2D electron motion plane is created 
only in a finite region — in a circle of radius r$. We 
showed that such a magnetic field fails to confine elec- 
trons both for the standard parabolic dispersion law or 
the ultra relativistic linear dispersion law for Dirac elec- 
trons in graphene. Although in such a magnetic dot no 
confined states are found, quasi-bound sates with a finite 
lifetime are present. 

An analysis of the quasi-bound states for the 
Schrodinger and Dirac electron was performed by means 
of the local density of states, and the position and 
width of the resonance peaks (the analogs of the quasi- 
stationary states), on the dot radius (or the strength of 
the magnetic field) were calculated. 

The broadening of these peaks (the inverse lifetime of 
the quasi-bound state) is mainly caused by the touching 
of the quantum dot border by the electron wave function 
tail. Due to the exponential character of this tail there 
exist a rather sharp border between broadened quasi- 
bound states and those which can be considered as nearly 
bound ones. This border was found by applying the com- 
plex energy eigenvalue method which is shown to be in 
qualitative agreement with the results obtained from the 
local density of states technique. 

It is shown that the difference of the quasi-bound states 
in the magnetic dot between the Schrodinger and Dirac 
electrons is only in the energies of these states which is a 
consequence of the different energies of the corresponding 
Landau levels. 

There is a single exception: in the case of Dirac elec- 
trons there exists a zero energy bound state for negative 
values of the angular momentum (the momentum which 
is opposite to the direction of the classical electron ro- 
tation along the Larmor circle). When the dot radius 
ro (or the magnetic field strength) decreases the degen- 
eracy of this zero energy level decreases skippingly while 
all other quasi-bound states disappear smoothly via their 
broadening. 
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